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Abstract 

We study the ion response to an initial perturbation in a weakly ionized plasma with ion flow 
driven by a dc electric field. The analysis is made by extending the classical Landau work [J. Phys. 
(USSR) 10, 25 (1946)] to the ion kinetic equation including ion-neutral collisions and a dc electric 
field. We show, in particular, that the complex frequencies of ion waves can be directly found from 
a known expression for the ion susceptibility [Phys. Rev. E 71, 016405 (2005); Phys. Rep. 27, 
997 (2001)]; this is not obvious from its original derivation, because it only aims to describe the 
ion response for real frequencies. 

PACS numbers: 52.30.-q, 52.25.Ya, 52.25.Mq, 52.25.Dg, 52.35.Fp, 52.35.Qz, 52.27.Lw 
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I. INTRODUCTION 



Low-pressure gas discharges are characterized by the presence of strong ion flow driven 
by the electric field that naturally arises to maintain the balance of absorption of ions and 
absorption of electrons on the electrodes and walls of the discharge chamber UQ|. This field 
extends far beyond the pure ion sheath , so that the velocity of the ion flow driven by 
this field often exceeds the thermal velocity of neutrals in the most part of the discharge 
(see, e.g., Fig. 7 of Ref. |4j]). 

This field-driven ion flow plays an important role in variousphenomena in low-pressure 
gas discharges, with the most obvious examples being ion waves [z-9] and interaction between 



charged dust particles levitated in discharges 



!CH17j. To study the role of the ion flow in 



such phenomena, it is helpful to derive the ion susceptibility in a plasma in the presence of 
field-driven ion flow. 

The ion susceptibility in the presence of field-driven ion flow has been self-consistently 



derived in Refs. 



18 



191 ] . To clarify, the self-consistency here means two things. First, the 



steady-state velocity distribution is found from the model itself, i.e., from the balance of ion 
acceleration in the dc electric field driving the ion flow and ion-neutral collisions (instead of 
assuming a model distribution, e.g., a shifted Maxwellian distribution). Second, collisions 
and the dc field not only define the steady state but are also accounted for in the analysis 
of perturbations. The resulting susceptibility is conveniently expressed in Ref. [l8j via the 
plasma dispersion function [E g. (8) of Ref. [18j]. 



The derivation of Refs. 



18 



19j aims to describe the ion response for real frequencies and 



real wave numbers and does not allow using the derived susceptibility to find the complex 
frequencies of ion waves (for real wave numbers), and the present paper addresses this issue. 

19[ does not allow using the derived 



18 



Let us first explain why the derivation of Refs. 
susceptibility to find the complex frequencies of ion waves 
We first note that the derivation of Refs 



18 



191 ] assumes the perturbations to be oc 



exp(— iut + ik • r), substitutes this to the ion kinetic equation, and finds the susceptibility 
as the normalized ratio of the complex amplitudes of the ion density and the potential. [Of 
course, for real u and k the perturbations of the above form can be considered as the Fourier 
components of the actual perturbations (as the Fourier transform is defined for real u and 
k), so that for real u and k the susceptibility derived by the above method has a clear 
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meaning: it is the normalized proportionality coefficient between the Fourier components of 
the perturbations of the ion density and the potential.] 

We point out the following well-known fact: in a collisionless one-component Maxwellian 
plasma in the absence of extraneous charges there are no solutions where the perturbation of 



the distribution function and the pertur 
the same complex u and the same real k 



Dation of the potential are oc exp(— iut + ik ■ r) with 



201 ] . (A one-component plasma is the approximation 



where only one plasma component oscillates.) Indeed, existing solutions in a collisionless 
one-component Maxwellian plasma in the absence of extraneous charges are characterized 
by different time dependencies of the potential and the distribution function (for the same 
k): while oscillations of the potential are exponentially Landau damped, the perturbation of 
the distribution function experiences undamped oscillations due to free dispersal of particles 
(oc exp[— i(k-v)i] at large t, see Ref. |20|, §1, last paragraph). Such solutions cannot be found 
when the perturbation of the potential and the perturbation of the distribution function are 
assumed to be oc expi—icut + ik ■ r) with the same uj and the same k (as assumed in the 



derivation of Refs. 



18 



. Il9j). Thus, it is not clear how the susceptibility of Refs. [lg, Il9j ) can 
be used to characterize such solutions. 

Leaving aside the issue with the assumed form of the perturbations, we also point out that 



the validity of the mathematical calculations of the derivation of Refs. [18|, [19| is limited to 
the case where the perturbations do not decay faster than the ion-neutral collision frequency 
[i.e., where \m.(uj) is larger than minus the ion-neutral collision frequency], as shown in 
Appendix [A] 



To find the complex frequencies of ion waves in the mode 
gest analytically continuing the susceptibility of Refs. 



18 



of Refs. 



19[ , one might sug- 



191 ] to the lower half of the complex 



frequency plane, in analogy to a collisionless Maxwellian plasma [20], but this requires sub- 
stantiation. Indeed, the correctness of this procedure for a collisionless Maxwellian plasma 



has been justified by solving the initial value problem 20 1, and it is not known how the so 



lution of the initial value problem for the model of Refs. 



18 



is not clear whether the same method (as in Ref. |20J|) of calculation of the integral restoring 



191 ] looks like. In particular, it 



18 



J9k 



18 



19]. 



the potential from its Laplace transform in time will work for the model of Refs. 

To address the issue, we solve the initial value problem for the model of Refs. 
That is, we consider an arbitrary initial perturbation and analyze its time evolution. This 
allows us to confirm the above hypothesis that the complex frequencies of ion waves can be 
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found by analytically continuing the susceptibility of Refs. [18|, |l9| to the lower half of the 



complex frequency plane. Since in Refs. 



18 



19| the derived susceptibility is presented in 



such a form that its analytical continuation to the lower half of the complex frequency plane 
is given by the same expression, it can be said that the complex frequencies of ion waves 
can be found from the same expression for the susceptibility as presented in Refs. 
by simply considering this expression as a function of complex u. 

19] has already been used in Ref. [2l| to study ion waves, 



The susceptibility of Refs. 



18 



and thus our work substantiates the validity of the study reported in Ref. 21(. We further 
note that, as our derivation can be easily generalized to a multiple- ion-species plasma, our 
work is of importance to the potential use of the susceptibility of Refs. [lg, Q to study ion 
waves in the presence of multiple ion species. 



II. MODEL 



We consider a weakly ionized plasma in a dc electric field E driving ion flow. We 
assume that electrons obey a Boltzmann distribution with a very large temperature, so 
we consider their number density uq to be homogeneous and fixed (note that Boltzmann 
electron distributions, in the presence of field-driven ion flow, are common in low-pressure 
gas discharges {§, 9|; we address the role of the electron temperature in Sec. [V}. We assume 
E to be homogeneous and use the kinetic equation for ions with the Bhatnagar-Gross-Krook 



(BGK) ion-neutral collision term 

Of 



18 



22 



dt dr m 



—uf + z/$ 



23[ and Poisson's equation: 
E 



df_ 

dr J dv 

/(v')rfv', 



-Ac 



e 



fdv-n 



(1) 



(2) 



where / is the ion velocity distribution function, <p is the electric potential describing the 
time-space varying field (i.e., the field apart from E ), 



exp 



2v 2 



(3) 



M (2vr<) 3 / 2 

is the normalized Maxwellian velocity distribution of neutrals, v is the ion-neutral collision 



frequency, which is assumed to be velocity-independent, v tn = ^Jh^Ljui is the thermal 



velocity of neutrals, T n is the temperature of neutrals, e is the elementary charge (ions are 
assumed to be singly ionized), m is the ion mass, k-g is the Boltzmann constant, and eo is 
the permittivity of free space. Note that the BGK term exactly describes charge transfer 
collisions under the assumption of a velocity-independent collision frequency, as explained 
in Ref. [18 1. 



The homogeneous steady-state solution / = f is found from Eqs. ([I]) and ([2} by setting 
= 0, df/dt = 0, df/dr = 0. The resulting velocity distribution f is not a shifted 



Maxwellian distribution and can be written as 



id 



where 



v f = ^ (5) 
mv 

and the subscript "f" stands for "flow", as the flow velocity (1/no) / v/odv can be shown 
to be equal to v f . Equation @ here is Eq. (3) of Ref. [18j, but rewritten using another 
integration variable in order to show that fo is an integral superposition of shifted Maxwellian 
distributions with exponential weights. 

We consider a small arbitrary initial perturbation of the ion distribution function and 
determine the resulting evolution of the potential cf). It is sufficient to consider the pertur- 
bations to be of the form 

/(r, v, t) - / (v) = /i(v, t) exp(zk ■ r), 

0(r,t) = <t>i(t) exp(ik • r), (6) 

(with a real k), since in the general case one can expand the perturbations into a Fourier 
integral in space and then consider each Fourier component separately. We denote /i(v, 0) 
by /i,i(v) (where the subscript "i" stands for "initial") and determine 4>i(t). 



III. RESULT 



In this section we formulate the result, while its derivation is provided in the next section. 
We first note that there are two cases, namely whether the initial perturbation satisfies 
certain reasonable conditions (stated in Sec. HVl) or not. These conditions are satisfied 
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for a wide class of functions /i,i(v); an example of an initial perturbation satisfying these 
conditions is /i,i(v) oc exp[— v 2 /(2w 2 )], where w is an arbitrary real constant. 

If the initial perturbation satisfies these conditions, then the evolution of the potential is 
described by: 



x (t) = C j exp(-iw i t) + H(t, T) exp(lY), 

Im(ajj)>r 



(7) 



where Cj are certain coefficients [for their explicit form, see Eq. (|32|) ]. the complex frequencies 
Uj are the roots of the dispersion relation 



o. 



(8) 



the index j numbers these roots (there is an infinite number of these roots, even when E = 
and v — > 0, and these roots are known as the Landau plasma modes 24J|). the sum is only 
over those j that satisfy lm(ojj) > T, F is an arbitrary real number (e.g., it can be chosen 
to be highly negative so that the exp(rt) factor in the second term of the right-hand side of 
Eq. (jTj) decays with time very rapidl y), the function x( w ) is given by the same expression 



as the susceptibility derived in Refs. [18 



19j: 



uj 



pi 



B{u) 



A(u) 
B(u) = 



v 2 1-A(w)' 

exp[— ^(u, rj)) dr], 
o 

00 ?7exp[-^(cu,77)] 



1 + z(k • v f /V)?7 
ioo 



dr), 



(9a) 
(9b) 
(9c) 



1 

+ 2 



ik ■ Vf 



+ 



v 



1 - 

V 



(9d) 



H(t, T) is a certain function [for its explicit form, see Eq. fl33|) ] and tends to zero as t — > oo 
(as shown in the last but one paragraph of Sec. HVl) . i.e., at large t the second term of 
the right-hand side of Eq. (J7]) is negligible as compared to any non-zero term of the sum 
represented by the first term of the right-hand side of Eq. (J7|, and to pi = [n e 2 /(e m)] 1 / 2 
is the ion plasma frequency. Equation (j7j) shows that the asymptotic behavior of <f>i(t) at 
large t is independent of the properties of the initial perturbation, as the frequencies Uj are 
determined entirely by the plasma parameters and the wave number k. The frequencies Uj 



are analyzed in Ref. 



2l|. 



If the initial perturbation does not satisfy the conditions mentioned in the beginning of 
this section, then the asymptotic behavior of <pi(t) at large t may be determined not only 
by the plasma parameters and the wave number but also by the properties of the initial 
perturbation, in contrast to the case considered above. 

Note that we also provide an expression for the potential </>i(£) valid irrespectively of 
which of the above two cases holds; this expression is Eq. (f3Tj) . 



IV. DERIVATION 



To derive the result stated above, we need to linearize Eqs. flTJ) and (J2]) with respect to f\ 
and 0i and then solve the resulting linearized equations with respect to <\>\. The linearization 
gives: 

dfi eE df\ ikefa df 

— + ik ■ v/i + -= -t— = 

at m av m ov 



-ufr + uQu j /i(v')rfv', (10) 



^ 2 01 = -//lrfv. (11) 

eo J 

We solve Eqs. ( jTOj) and ( TTTT) by using the Laplace transform in time. The Laplace trans- 
form in time is defined here as 

poo 

$(u)= / ip{t) exp(iwt) dt, (12) 
Jo 

where the hat denotes the application of the transform and if)(t) is an arbitrary function. 
The inverse transform restoring ip(t) from ip(uj) is: 

i /•+oo+t7 
-oo+«7 

where the integration is performed in the complex plane over the line Im(u;) = 7, and 7 is any 
real number such that this line is in the region of convergence of the integral in Eq. (1T21) . 
Taking the Laplace transform of Eqs. ffTD]) and (TTTT) . i.e., multiplying these equations by 
exp(icot) and then integrating over J °° dt, we get 

(-zu; + v + ik ■ v)/i H — = — 

m ov m av 



1 r+oo+i-y 

ip(t) = — / t/j(u) exp(-iut) du, (13) 

^ J -oo+i'y 



+z^m J fiWdv' + fo, (14) 



k 2 ^ = - Adv. (15) 
£o J 

Thus the initial value problem is reduced to solving Eqs. (|I4]) and ( JTBT) with respect to 0i. 

Let us reduce solving Eqs. ( IT4l) and ( JT5l) with respect to 0i to finding the solution of 
Eq. f lT^|) with respect to /i for a fixed </>i. To do so, we first note that the latter solution 
has the following structure: 

A = ^101 + L 2 J A(v') dv' + C(fo), (16) 

where the first, second and third terms in the right-hand side denote the direct contributions 
from the first, second and third terms of the right-hand side of Eq. ffT^j) . respectively, Li j2 are 
certain coefficients (which depend on v), £ is a certain linear operator, and we do not need 
explicit expressions for L 12 and £ at the moment. [Note that the solution ffTB"]) is unique 
because of the boundary condition fi = 0.] We integrate Eq. ( JT6l) over v, which allows 
us to express the integral f fi dv via Li j2 , £, and <f>i. By substituting the result to Eq. ( TT5l) 
and using Eq. ( TT3j) . we get 

-+00+17 



0i (t) = 7^- du exp 



-icut) 



(e/e ) f £(fi,j)dv 
k 2 [l-jL 2 dv]-(e/e )fLidv' 1 J 

Thus, to solve the initial value problem, there only remains to derive explicit expressions for 
L 12 and £. 

Let us now solve Eq. f|LT|) with respect to f\ for a fixed 0i in order to derive the required 
expressions for Li j2 and £. Making the substitution 

A=/3exp(-a), (18) 

where 



m 

a 



{-iu + ik x v x + v)v z + -i^u^ 



(19) 



eE 

the z-axis is in the direction of E , and the x-axis is directed so that k y = 0, we reduce 
Eq. (dU) to 

d(3 _ mexp(q) [ 

where [. . .] denotes the right-hand side of Eq. (|T^|) . To integrate Eq. (|20|) . we need a boundary 
condition on 0, so let us find a boundary condition on /3 from that fi = at v z — > ±00. Let 
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us explicitly write the condition that j\ = at v z — > ±00 in terms of /3: /3exp(—a) = at 
v z — > ±00. The behavior of exp(— a) at v z — > ±00 is determined by the sign of Im(o;) + v 
[as seen from Eq. (TI9"]) ]. but we only need to consider the case 

Im(w) + v > (21) 

because in Eq. (j!7p we can choose any 7 larger than a certain threshold. Considering this 
case, we note that exp(— a) does not tend to zero as v z — > —00 and tends to zero as v z — > +00, 
as seen from Eq. f fT9|) . Thus, we obtain that /3 = at v z — > —00; note that one cannot set 
(3 = at v z — > +00. Using that j3 = at v z — > —00, we integrate Eq. (120]) to obtain 



Li — — I dv' z exp(a' — a 



f - ul) k z — - -^-Jo 

1 



(22) 



L 2 = — / dv' exp(a' - aW u , (23) 



00 

Vz 



m 

C(fu) = — I dv' z exp(«' - a)f 1A , (24) 

where all quantities denoted by the prime symbol contain v' z instead of v z . 

Thus we have in principle solved the initial value problem, as the solution is provided by 
Eqs. ( TlTI) and (I22"l)-(l24"l). with 7 being a sufficiently large number such that the result of the 
integration in Eq. (JTTJ) for any larger 7 is the same. 

Let us obtain a simpler form for the potential by simplifying the integrals j L 12 d\ and 
J £(fi ti ) dv present in Eq. ffTT|) . 

Let us start with the integrals j L± t 2 dv. We first perform the integration over v y and 

then over v x using the formula 

+00 / 2 \ / „2 

2.1,2 „„„ / X , \ J„ _ a„1,2 r— I 

00 



x ' exp ( — — + cx j dx = ic ' v2tt exp ( — — ) (25) 



valid for complex c. The next step is to replace v' z by v z — VfT), where rj is a new inte- 
gration variable, and to integrate over v z using Eq. ( )25l) . Then the integration over £, 
the variable originating from Eq. (jlj), is elementary. These calculations yield j L\dv = 
— (k 2 n eB)/(mi> 2 ) and J L 2 dv = A, where A and B are defined by Eqs. f l9b[) and f )9c|) . 
respectively. Thus, the denominator in Eq. ( fl7|) can be written as (1 + x)/c 2 , where x is 
given by Eq. ©. 
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Let us now simplify the integral J £{fi,i) dv. We first change the order of integration, 
dv z dv' z —> dv' z dv z . Then we replace v z by v' z + q, where q is a new 
integration variable, and integrate over q using the formula 

i — 2 sign (6) nr 

2 



exp ( —ax — -ibx 2 I dx 



x exp 



[i sign(6) — l]a 



x 



1 + erf 



[i sign(fe) — l]a 



(26) 



valid for real a > and real b ^ 0, where 



erf(*) = 4= f exp(-£ 2 )^ 
V 71 " Jo 



(27) 



is the error function defined for a complex variable. After that we can remove the prime 
symbol in v' z . The result can be written as: 

sign(fc z ) + i 



£(/i,i)dv 



m 



eE \k z 



h,iZ{K) dv, 



(28) 



where 



Z(x) = 2i exp(— x 2 

is the plasma dispersion function and 

sign(& z ) + i 



exp(-£ 2 )d£ 



(29) 



A 



rn 



2 V eEo\k z \ 
x (— iu + v + ik • v). 



(30) 



Using the above obtained results for J Li )2 rfv and J £(/i,0 dv, we write the potential 
13) as 



01 (*) 

+00+17 



sign(A; z ) + i / me 



4ne k 2 \J E Q \k z \ 
j / 14 Z(A) rfv 



x / gL> exp(— zu;i) 



(31) 



where x is given by Eq. OH]); again, we note that 7 is a sufficiently large number such that 
the result of the integration in Eq. ( 131 j) for any larger 7 is the same. 



10 



Im(c£>) 



Re(ffi) 



























FIG. 1: Contour of integration used to derive Eq. (JT]). The dots denote the solutions of the 
dispersion relation (jSJ). 



Let us now calculate the integral in Eq. (1311 by using Cauchy's integral theorem and 
assuming certain reasonable conditions on the initial perturbation. To do so, we consider 
the fraction in the last line of Eq. fl3T|) as a function of u in the entire complex u plane 
[despite that we only obtained the expressions for L 12 and C for the region (121 p ] and modify 
the contour of integration over u in such a way that (i) Cauchy's integral theorem guarantees 
that the value of the integral is not changed and (ii) integration over the new contour yields 
Eq. ©. 

We assume that the initial perturbations satisfies the following conditions: (i) the ex- 
pression in the numerator in the last line of Eq. (13T|) defines an analytic function of u in the 
entire complex u plane and (ii) this function is sufficiently well-behaved at Re (a;) — > ±oo 
so that there are no problems with the validity of the results stated below. We verified 
numerically that conditions (i) and (ii) above are met for a wide class of functions /i,i(v) 
[e.g., for /i,i(v) oc exp(— v 2 /2w 2 ), where w is an arbitrary real constant]. We displace the 
integration contour down as shown in Fig. [TJ This, according to Cauchy's integral theorem, 
does not change the value of the integral, as the function 1/(1 + %) [where x is defined by 
Eq. (Q in the entire complex u> plane] is analytic in the entire complex u> plane except at the 
zeros of the function 1 + X- The integral over the new contour is the sum of the contribution 
of the circles shown in Fig. [Hand the contribution of the horizontal part. This allows us to 
write Eq. ([7]), where 

. — i sign(fc 2 ) / me J fi^Z^A) dv 



2e k 2 



E \k z 



dx/du 



(32) 



11 



H(t,T) 



sign(&; 2 ) + % I me 



x 



— oo 



oo 



dw r exp 




(33) 



and the number T defines the position of the horizontal part of the new contour as shown 
in Fig. [TJ we emphasize that in Eqs. fl32|) and f )33|) . x an d A are given by Eqs. @ and f )30|) 
m t/ie entire complex u plane despite that we only obtained the expressions for L 12 and 
C for the region ( 1211) ; the first term of the right-hand side of Eq. (I7j) is the contribution 
of the circles shown in Fig. [TJ and the second term of the right-hand side of Eq. ([7]) is the 
contribution of the horizontal part; we introduced a new integration variable u r = Re(o;) 
(where the subscript "r" stands for "real") to conveniently write the contribution of the 
horizontal part; the integration over u r in Eq. (|33|) is performed over the real axis. 

To substantiate our claim made in Sec. IIHI that the function H (t, T) tends to zero as 
t — > oo, we note that in the right-hand side of Eq. ( 1331) . t is only present in exp(— icj r t), 
so that the integral in Eq. ( )33|) can be considered as the inverse Fourier transform in time. 
Hence, this integral tends to zero as t — > oo, according to the Riemann-Lebesgue lemma 
[indeed, the numerator in this integral is a sufficiently well-behaved function, as assumed 
in the beginning of the previous paragraph; the denominator in this integral tends to 1 as 
u T — > ±oo, as can be easily verified]. 

Let us now briefly discuss the case where the initial perturbation does not satisfy the 
conditions stated in the beginning of the second paragraph back. To deal with this case, one 
can analytically continue the numerator in the last line of Eq. ( 13~T|) into the region below the 
line Im(o;) = 7. To do so, one may need to specify cuts; the obtained analytical continuation 
may also have singular points. Then one can attempt to modify the procedure detailed 
in the second paragraph back by modifying the contour shown in Fig. [TJ in particular, 
the contour should now also pass around the aforementioned cuts and/or singularities so 
that Cauchy's integral theorem can be applied. The contributions from the contour parts 
that pass around the aforementioned cuts and/or singularities should be added to Eq. (IT]). 
The time dependence of these additional terms will be determined not only by the plasma 
parameters but also by the properties of the initial perturbation. 
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V. ROLE OF THE ELECTRON TEMPERATURE 



Let us discuss the role of the electron temperature. There are two effects related to a 
finite electron temperature: (i) the electron response to ion oscillations, and (ii) a finite in- 
homogeneity length of the Boltzmann electron distribution in the field E . Concerning effect 
(i), it can be taken into account by adding the term l/{k\ e ) 2 to the left-hand side of Eq. (jSj), 
where A e = [e k B T e / (n e 2 )] 1//2 is the electron Debye length and T e is the electron temper- 
ature. As regards effect (ii), the corresponding inhomogeneity length is L e = k-QT e / (eE ), 
and our model applies when this distance is larger than both the ion-neutral collision length, 
which is maxjuf, v tn }/u, and the wavelength 2ir/k. This means the following applicability 
limitation of the model: 

— — > max jv f , VtgVt, — | . (34) 



VI. CONCLUSION 



We have demonstrated that the expression for the ion susceptibility derived in Refs 



.0 



19] [and given by Eq. (jUj) of our paper] can be used in the entire complex frequency plane to 



find the complex frequencies of ion waves, despite that its original derivation 
aims to describe the ion response for real frequencies. 



18 



19( | only 
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Appendix A: Note on the derivation of Refs. 



HQ 



In this Appendix we substantiate our claim made in Introduction that the validity of the 
mathematical calculations of the derivation of Refs. 
0. To do so, we repeat the derivation of Refs. 



18 



18 

3 



19] is limited to the case lm(u) + u > 



and point out where this restriction 
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arises. In this Appendix we use the notations of the present paper and not the notations of 



18 



19] are the same. 



Refs. [La Ll9fl: the model of our paper and the model of Refs. 

nn 

The idea of the derivation of Refs. [18|, |19J is to assume that the perturbations are of the 
form 

/(r, v, t) - /o(v) = A(v) exp(-iujt + ik • r), (Al) 
</>(r, t) = a exp(— uot + ik ■ r), (A2) 

(here the subscript "a" stands for "amplitude"), substitute this to the kinetic equation 
linearized with respect to the perturbations, solve this equation with respect to / a for a 
fixed a , and find the susceptibility as 

x= <U£l. (A3) 

tok 2 (j) a 

Concerning the first step, the linearized kinetic equation takes the form: 

w eE df a eik(p a df 

+ v + ik ■ v)/a H — = — 

m ov m av 

[ f a (y')dV. (A4) 



To solve the linearized kinetic equation [Eq. (1A4j) ] with respect to / a , one makes the 
substitution 

/ a = i 71 exp(-a) (A5) 

where a is defined by Eq. ( fl9|) . The resulting equation is 

dF _ mexp(q) 

^"^^ [ --- J ' (A6) 
where [. . .] denotes the right-hand side of Eq. ( 1A4I) . 

To integrate the resulting equation [Eq. ( 1A6I) ]. one needs a boundary condition on F, 
and here there are two different cases, namely whether Im(o;) + v is positive or negative, as 
explained in the following. The boundary condition on / a is 

/ a = at v z ±oo. (A7) 

Let us explicitly write this condition in terms of F using Eq. (1A5I) : 

Fexp(-a) =0 at v z — »■ ±oo. (A8) 

The behavior of exp(— a) at f z — >■ ±oo depends on the sign of Im(co>) + u, as seen from the 
definition of a [Eq. ( fl9]) ]. so there are two different cases: 
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1. lm(u) + v > 0. In this case exp(— a) does not tend to zero as v z — > — oo and tends 
to zero as v z — >• +00, as seen from Eq. f|T9|) . Using this fact and Eq. flA8j) . we get the 
following boundary condition on F: 



at v 7 — t- —00 



(A9) 



(i.e., one cannot set F = for t> z — > +00). Having obtained Eq. (IA9[) . we do not need 
to know the behavior of F at — > +00 in order to integrate Eq. ( ]A6j) . as Eq. ( 1A6I) is 
a first-order differential equation. 

2. Im(o;) + z/ < 0. In this case, in contrast to the preceding case, exp(— a) tends to zero as 
v z — > —00 and not as v z — > +00, as seen from Eq. ( IT9l) . Thus we obtain the following 
boundary condition on F: 

F = at v z -+ +00 (A10) 
(i.e., one cannot set F = for f-, — >• —00). 



The derivation of Refs. 
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19[ only considers case 1 above and is thus limited to the case 



Im(a;) + v > 0. Using Eq. flA9j) . one integrates Eq. flA6j) to obtain: 



fa 



dv z exp(a' — a) 



v xk x f/\, 



/a(v") rfv" 



(All) 



where all quantities denoted by the prime symbol contain v' z instead of v z , and the x-axis is 
directed so that k y = 0. Equation (IA11I) above is Eq. (5) of Ref. [19J [one needs to clarify 
here that Eq. (5) of Ref. 191 is incorrectly typed, as the quantities inside the large round 
brackets in Eq. (5) of Ref. [19] should contain v' z instead of v z ] . 

Let us demonstrate that Eq. flAHI) does not apply in case 2 above. To demonstrate this, 
in the following we show that in case 2 Eq. (lAllj) does not satisfy the boundary condition 
/ a = at v z — > +00. We consider the particular case k z = (for simplicity) and rewrite 
Eq. (lAlip so as to explicitly show the dependence of / a on v z at v x = 0: 



/ a ex exp 
dv' z exp 



mv z 



-iuj + v) 



mv' 



,'2 



eEn 



—ioj + v) — 



2v 2 



(A12) 
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As v z — > +00, neither the integral in Eq. (1A12j) nor the exponent in front of this integral 
tends to zero, so we see that the condition f a = at v z — > +00 is indeed violated. 
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